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Abstract 


We report on long-term multi-wavelength monitoring of blazar Mrk 421 from 2010 to 2011. The source 
exhibited extreme X-ray flares in 2010. Our research group performed optical photopolarimetric follow-up 
observations using the Kanata telescope. In 2010, the variability in the X-ray band was significant, while 
the optical and ultraviolet (UV) flux decreased gradually. Polarization properties also exhibited unique 
variability in 2010, suggesting the presence of systematic component of polarization and magnetic field 
alignment for the emergence of a new polarized emission region. In contrast, in 2011 the variability in the 
X-ray band was smaller, and the variability in the optical and UV bands was larger, than in 2010. To 
explore the reasons for these differences, spectral fitting analysis was performed via simple synchrotron-self 
Compton modelling; the results revealed different behaviors in terms of spectral evolution between these 
periods, suggesting different variability mechanisms between 2010 and 2011. In 2010, the radiation was 
likely the result of energy injection into the emitting regions with an aligned magnetic field. In contrast, 
in 2011 the superposition of different emission regions may have contributed to the low degree of observed 
polarization. It also implies that high-energy electron which were not accelerated to ultra-relativistic 
velocities were injected in 2011. 


1. Introduction 

Blazars are highly variable active galactic nuclei (AGN) 
that can be detected for all wavelengths, ranging from ra¬ 
dio waves to gamma rays. The strong relativistic jets from 
blazars are aligned with the observer’s line of sight and ap¬ 
pear bright due to relativistic beaming, having a luminos¬ 
ity typically greater than 10^^ erg s“^ (Urry & Padovani, 
1995). Blazar emission typically consists of two spectral 
components: one is attributed to synchrotron radiation 
from relativistic electrons, emitted at lower energies peak¬ 
ing in the radio through optical regime; the second occurs 
at higher energies from inverse Compton scattering radia¬ 
tion, peaking in the gamma-ray regime. A notable charac¬ 
teristic of blazars is their rapid, high-amplitude intensity 
variations or flares. Blazars can be classified into three 
types, based on the peak synchrotron radiation frequency 
i^peak (Ackermann et ah, 2011): low-synchrotron-peaked 
blazars (LSP), intermediate-synchrotron-peaked blazars 
(ISP), and high-synchrotron-peaked blazars (HSP). Due 
to relativistic effects, the radiation from jets dominates 
the overall spectral energy distribution. As a result, the 
optical band spectra of blazars are essentially featureless, 
compared with other AGNs. This makes blazars highly 
suitable for jet study. 

Polarized radiation from blazars, an indicator of low- 
energy synchrotron radiation, can vary sharply. Thus, 
the polarization of blazars is of interest for understanding 


the origin, confinement, and propagation of jets, due to 
the dependence of the polarization on the jet’s magnetic 
field structure (e.g., Visvanathan & Wills, 1998). Mead 
et al. (1990) performed a large-sample study of blazars in 
the optical band and found that a high polarization de¬ 
gree (PD) and variability of polarization are common phe¬ 
nomena in blazars. Ikejiri et al. (2011) reported statisti¬ 
cal photopolarimetric observations of blazars with a daily 
timescale, and suggested that the lower luminosity and 
higher peak frequency of synchrotron radiation objects 
(e.g., an HSP blazar) have smaller amplitudes for varia¬ 
tions in the flux, color, and polarization degree. However, 
questions remain regarding the reason for the difference 
between HSP and LSP blazars. What seems to be lacking 
is simultaneous multi-wavelength and polarimetric obser¬ 
vations. Multi-wavelength observation is valuable because 
the time-resolved spectral energy distribution (SED) in 
the order of physically relevant timescales allows verifi¬ 
cation of a synchrotron-self-Gompton (SSG) model that 
considers the time-development of the electron spectra. 
Wide-range, temporally deep, multi-wavelength polariza¬ 
tion studies should reveal some of the blazar’s character¬ 
istics, such as its energy dependence, time variability, and 
jet structure. 

Markarian 421 (Mrk 421; R.A. = 11^04™27.3^ Deck 
= -b38° 12' 31".8, J2000, z=0.031), one of the most 
observed blazars, is classified as an HSP blazar (Abdo 
et al., 2010b). Additionally, it was the first extragalactic 
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source detected at TeV energies, over the range 0.5-1.5 
TeV, by the Whipple telescopes (e.g., Punch et ah, 1992). 
The SED exhibits a double-peaked feature with two max¬ 
ima; its variations are basically well-understood within 
the framework of a simple one-zone SSC model (Band & 
Grindlay, 1985; Abdo et ah, 2011). Hence, Mrk 421 is a 
good target for flare mechanism studies about relativistic 
jets. We monitored extreme activity in the Mrk 421 X-ray 
bands that occurred in 2010, via long-term photopolari- 
metric measurements with the Kanata telescope. Here, 
we present our findings about the relationship between 
the variability observed in the synchrotron spectra and 
optical polarization measurements during the 2010 X-ray 
flare. We assumed a flat ACDM cosmology, with Hq = 75 
km s“^Mpc“^ and = 0.3 (Wright, 2006). 

2. Observation 

2.1. Observation using the Kanata telescope 

We performed V and i?c-band photometry and po- 
larimetry observations of Mrk 421 from December 2009 to 
May 2011, using the HOWPol instrument installed on the 
1.5-m Kanata telescope located at the Higashi-Hiroshima 
Observatory, Japan (Kawabataet ah, 2008). We obtained 
74 daily photometry data points for the V band, and 
68 photometry and polarimetry data points for the Rc 
band. A sequence of photopolarimetric observations con¬ 
sisted of successive exposures at four position angles of a 
half-wave plate: 0°,45°,22.5°,and 67.5°. The data were 
reduced under standard charged coupled device (CCD) 
photometry procedures. We performed aperture photom¬ 
etry using APPHOT packaged in PYRAF and differential 
photometry with a comparison star imaged in the same 
frame as Mrk 421. The comparison star is located at R.A. 
= 11:04:18.2 and Deck = -1-38:16:30.5 (J2000); its mag¬ 
nitudes are V = 15.57 and Rc = 15.20 (Villata et ah, 
1998). The data were corrected for Galactic extinctions 
A(V) = 0.050 and A(R) = 0.040. 

Note that it is also important to correct the instru¬ 
mental polarization caused by reflection at the tertiary 
mirror of the telescope. The instrumental polarization 
was modeled as a function of the declination of the object 
and the hour angle of observation, which was subtracted 
from the observed value. We confirmed that the accuracy 
of the instrumental polarization subtraction was smaller 
than 0.5% for the Rc band, using unpolarized standard 
stars. The polarization angle (PA) was defined in the stan¬ 
dard manner (measured from north to east), based on cal¬ 
ibrations with polarized stars, HD183143 and HD204827 
(Schulz & Lenzen, 1983). Because the PA has an ambi¬ 
guity of ±180° X n (where n is an integer), we selected 
n, which gives the minimal angle difference from previous 
data, assuming that the PA changes smoothly. We also 
confirmed that the systematic error caused by instrument 
polarization was smaller than 2°, using the polarized stars 
as mentioned. 

^ PYRAF is a product of the Space Telescope Science 

Institute, which is operated by AURA for NASA 

http: //www. stsci.edu/institute/software_hardware/pyraf 


2.2. Observation using Swift UVOT and XRT 

Swift observed the source between 2009 December 15 
and 2011 February 22. A total of 125 data points in the 
UVW2 band, 122 data points in the UVM2 band, and 111 
data points in the UVWl band were obtained with the UV 
and Optical Telescope (UVOT). UVOT data were reduced 
following the standard procedure for CCD photometry. 
Counts were extracted from an aperture of 5 arcsec radius 
for all filters and converted to flux using the standard 
zero points. (Poole et ah, 2008). The data were corrected 
for Galactic extinctions A{UVW2) = 0.136, A(UVM2) = 
0.143, and A{UVW1) = 0.098 (Schlegel et ah, 1998). 

We also used data obtained by the X-Ray Telescope 
(XRT) in the same epoch. There were 165 data points ob¬ 
served by XRT. All XRT data presented here were taken 
in windowed timing (WT) mode^, due to the high flux 
rate of the source. Data reduction and calibration of the 
XRT data were performed with HEASoft v6.4 standard 
tools. Events in the energy range from 0.2 - 8 keV, with 
grades 0-2, were used for the analysis. Source spectra 
were binned to ensure a minimum of 20 counts per bin 
to facilitate minimization fitting. Response files were 
generated with xrtmkarf Task, with corrections applied 
for point-spread function (PSF) losses and CCD defects. 
Spectral analysis was performed with the XSPEC software 
package version 12. An estimate of the absorption at low 
energies, caused by interstellar gas, is required for spectral 
curvature analysis. Because there was no evidence of the 
host galaxy in Mrk 421 by optical observation (Urry et ah, 
2000), we used an absorption model having a hydrogen- 
equivalent column density (JSh) similar to that of our 
galaxy, Nh = 1.61 x 10^°cm~^ (Lockman & Savage, 1995). 
We performed spectral analysis by fixing the Nh absorb¬ 
ing column densities to the galactic values. In Massaro 
et al. (2004) and Tramacere et al. (2007), the spectral 
shape of Mrk 421 in the X-ray band was well fitted with 
a log-parabolic model, as opposed to a power-law form. 
Spectral fitting with the log-parabolic shape is described 
below: 

F{E) = ph cm"^ s"^ keV^ (1) 

where a is the photon index at 1 keV and b is the spectral 
curvature. Figure 1 shows the spectral fitting results for 
the power-law and log-parabolic models. A Xr f®st for 
the fit with a simple power-law model indicates a reduced 
value, Xr = 3.57548 (d.o.f = 687), for the data taken on 
MJD=55148. In contrast, the reduced Xr value using the 
log-parabolic model showed that Xr = 1-279 (d.o.f. = 686) 
for the same data. Improvement in the reduced Xr values 
was also evident for other observational data. Therefore, 
we used the log-parabolic model in this paper. 

2.3. Observation using MAXI 

We used archival data^ of Mrk 421 using Monitor of 
All-sky X-ray Image/Gas Slit Camera (MAXl/GSC) ob- 

^ see Burrows et al. (2005) for more details on XRT observing 

modes 

® http://maxi.riken.jp/top/ 
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Fig. 1. Swift X-ray Telescope (XRT) spectral analysis of Mrk 421 (MJD=55148). Left: Spectral fit using the power-law model 
with systematic deviations on both sides of the residuals. Right: Spectral fit with the log-parabolic model with no deviation. 


servations. GSC is designed to scan the entire sky every 
92-min orbital period in the 2-30 keV band, giving it 
the highest sensitivity of the all-sky X-ray monitors. The 
data used in this analysis were taken between December 
2009 and May 2011. The energy range was set to 2 - 20 
keV in 1-day data bins. Details of the data analysis of 
Mrk 421 for the 2010 January and February flares were 
reported in Isobe et al. (2010). 

3. Results 

3.1. Light Curves 

Figure 2 shows the multi-wavelength light curves from 
the optical to X-ray bands between MJD 55180 and 55630 
(corresponding to 2009-12-15 and 2011-03-10). Below we 
describe details of the light curves for each energy band. 

3.1.1. X-ray Flux and Photon Index 

Photon flux at 2 - 20 keV measured by MAXI/GSC 
and the energy flux at 0.3 - 8 keV measured by Swift/XKT 
showed similar trends in variability. In the first half of the 
observations (< MJD ~ 55400, hereafter referred to as ’pe¬ 
riod 1’), X-ray band variability was prominent compared 
with that in the last half (> MJD ^ 55400 ; ’period 2 ’). On 
MJD 55250, the X-ray flux reached 30.6^0^ x 10“^° [erg 
cm“^ s“^] at its maximum value. In period 2, however, 
the variability for the X-ray band was negligible. The av¬ 
erage flux in period 2 was F 2 ;ave = 1-5 x 10“^° [erg cm“^ 
s“^]. The maximum flux in period 1 and average X-ray 
flux values for period 2 differed by a factor of ~ 20. MAXI 
observations indicated that there were no large flares from 
MJD 55400 to 55530; however, the observation of Swift 
was limited by the location of Mrk 421, close to the sun. 
The photon index also showed clear variability in period 
1, indicative of the trend of ’harder when brighter’. A 
similar trend was reported in Tramacere et al. (2009) for 
this source for previous X-ray flares. Note that photon 
indices for period 2 were softer than the values for period 
1 . 


3.1.2. UV and Optical-band Flux and Color 

UV and optical fluxes also showed clear variability, but 
their trends were different from that exhibited by the X- 
ray flux. In period 1, flux variability gradually decreased, 
with no significant flares correlated with the X-ray flare. 
In period 2, the UV and optical flux changed significantly. 
The i^c-band maximum and minimum flux differed by a 
factor of ^ 3 in period 2, while the X-ray flux remained 
nearly constant. The UV and optical color were statisti¬ 
cally compatible and relatively constant. The spectrum 
shape variability of the UV and optical bands was smaller 
than that for the X-ray band. 

3.1.3. Optical Polarization 

With regard to polarization, the PD exhibited large 
variability in period 1. Similar to the variability of the 
X-ray flux, the variability in PD for period 1 was larger 
than that observed for period 2. In period 2, the PD was 
low and inactive. In contrast, the variability of the PA 
was small for period 1. The average value of PA for pe¬ 
riod 1 was ^ 150 degree. As opposed to the variability 
of PA observed in period 1, The PA in period 2 changed 
significantly. Such a rotation of the PA was similar to 
that observed for BL Lac and 3G 279 (Marscher et ah, 
2008; Abdo et al., 2010a); however, this observation re¬ 
quires closer examination, due to the pseudo-rotation of 
PA that appears when the PD is low. Figure 3 shows the 
polarization on the Stokes parameter QU plane for each 
period. Variation on the QU plane implies rotation in pe¬ 
riod 2; however, this variation could also be attributed to 
random motion about the origin of the QU plane. 

The variation on the QU plane suggests a difference 
in the polarization between periods 1 and 2. In figure 
3, the locations of Q and U were systematically differ¬ 
ent for periods 1 and 2: specifically, there was some 
systematic component of polarization in period 1. The 
average values of Q and U for each period were Q^ = 
0.013±0.005, Ui = -0.019±0.005, Q 2 = 0.000±0.005, and 
U 2 = —0.004 ± 0.005, respectively. To examine the vari¬ 
ability in more detail for period 1, we separated the data 
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Fig. 2. Light curves in multi-wavelength observation of Mrk 421 between Jan. 2010 and Feb 2011. The top 
panel shows 2-20 keV flux measured by MAXI/GSC. The second panel shows 0.3 - 8 keV flux measured 
by Swift/XRT. The third panel shows the photon index at 1 keV, measured by Swift/XRT. The fourth panel 
shows the UV and optical flux measured by Swift/WOT and Kanata/HOWPol. The fifth panel shows UV and 
optical color. The sixth panel shows the polarization degree in the Rc band and the PA (bottom panel). 


points in period 1 into three epochs, based on the X-ray 
light curve (also see figure 2): Pl-1 (MJD 55180-55227), 
Pl-2 (MJD 55227-55303), and Pl-3 (MJD 55303-55350). 
The right panel of figure 3 shows the detailed behavior of 
variability of QU in period 1. These results indicate that 
the variability of PA in each flare has a large variance. In 
other words, the PA did not appear to be completely con¬ 
stant during X-ray flares. We also checked for the presence 
of a systematic component of polarization, as reported by 
Morozova et al. (2014). for period 2, but no significant 
component was confirmed. 


4. Discussion 

We performed photopolarimetric observations of 
Mrk 421 in an X-ray active state between 2010 and 2011. 
Simultaneous wide-band spectra of synchrotron radiation 
were obtained on a daily timescale. Generally, HSP 
blazars, such as Mrk 421, tend to not have good corre¬ 
lation between the polarization degree and flux in the 
optical band (e.g., Ikejiri et ah, 2011). The correlation 
between the X-ray flux and optical flux is also unclear for 
HSP blazars (e.g., Finke et ah, 2008); sometimes they are 



















5 


No. ] 




Q [%] Q [%] 

Fig. 3. Polarization on Stokes QU plane for each period. In the left figure, red indicates period 1 (MJD 
55180-55430) and black indicates period 2 (MJD 55431-55630). In the right figure, the data points in period 1 
are separated into three epochs; (Pl-1, MJD 55180—55227; Pl-2, MJD 55227-55303; and Pl-3, MJD 55303-55350). 


correlated and sometimes not. Therefore, our results sug¬ 
gest that we should take into account the type of flare 
when investigating corresponding changes in the X-ray 
band flux. 

In the case that individual X-ray flares in period 1 
were produced by a single emission region, a constant PA 
and good correlation between the PD and flux were ex¬ 
pected. Our observations, however, indicated some varia¬ 
tion in the PA; additionally, no clear correlation was evi¬ 
dent between the flux and PD using a timescale of days. 
Systematic component of polarization was detected in pe¬ 
riod 1 (also see figure 3). Taken together, these results 
suggest that the individual X-ray flares may be the result 
of the superposition of a long-term component (timescale 
of weeks to months) and several short-term components 
(timescale of days), as opposed to a model using the super¬ 
position of short-term components without the long-term 
component (e.g., Moore et ah, 1982; Uemura et ah, 2010). 
This would explain the systematic difference of QU be¬ 
tween period 1 and period 2 , as well as the complexity of 
the correlation between flux and PD. 

The emergence of a new polarized emission region asso¬ 
ciated with the long-term component accompanied by an 
X-ray flare in period 1 also indicates that such extreme 
X-ray flares originate in a region where the magnetic field 
is aligned. Shocks in the jets lead to variations in the local 
emissivity and produce outbursts. Due to the compression 
of the magnetic fields during shocks, these variations will 
also lead to variations in the polarization. Additionally, 
the average PA of 150 degrees, without any large changes 
in PA such as the rotation or swing observed in period 2, 
is in good agreement with the position angle of the parsec- 
scale jet measured in the radio band (Piner & Edwards, 
2005). The conformity of the position angle of the jet and 
the PA have been shown for several blazars, suggesting a 
’shock-in-jet’ model (e.g. Hagen-Thorn et ah, 2008; Itoh 
et al., 2013). Namely, the magnetic field, compressed by 
shocks propagating down the jet, becomes aligned in a 
direction transverse to the jet. Note that long-term com¬ 
ponent of PA also suggests the presence of a persistent 


magnetic field in the underlying jet, because the timescale 
of the long-term quasi-stationary PA is somewhat long. 
Small variations in the polarization may be explained by 
the superposition of different subcomponents. Details of 
this scenario are discussed later in this section. 

In contrast to the behavior in period 1, the X-ray flux 
and polarization activity were low in period 2 , while the 
optical flux showed large variability. This implies that 
the variability mechanism was different between periods 
1 and 2. Thus, we attempted to resolve the difference 
between these periods using broad-band spectral fitting. 
The emission mechanism of X-ray flares is well studied by 
modeling the SED with the synchrotron and SSC model. 
Here, we researched the evolution of the SED for peri¬ 
ods 1 and 2 by spectral fitting with a simple SSC model. 
The simplest leptonic model typically used to describe the 
emission from an HSP blazar is the one-zone SSC model 
(Einke et ah, 2008). As mentioned above, the SED in pe¬ 
riod 1 may be consistent with a two-component model. 
However, for the sake of simplicity, we fitted the data us¬ 
ing a one-zone model for each period and compared the 
jet parameters of periods 1 and 2 . 

In the model, radio to X-ray emission is produced 
by synchrotron radiation from electrons in a homoge¬ 
neous, randomly oriented magnetic field (H), and gamma- 
rays are produced by inverse Compton scattering. We 
parametrized the electron distribution as a broken 

power law, expressed as 


iVe(7) = ^eX 



7 < 7& 
otherwise. 


( 2 ) 


where 7 is the Lorentz factor of the electron, 7 f, is the 
break energy of the electron, and Kg is a normalization 
term for the electron spectrum within the electron Lorentz 
factor from to ymax- Note that we also tried to fit 

the spectra using a single power-law model of the electron 
distribution, but the fitting results showed physically in¬ 
correct values (e.g., the Doppler factor was > 100 for the 
single power-law model). Similar results were reported in 
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Finke et al. (2008). Thus, we selected the broken power- 
law model for use in this study. Because our investiga¬ 
tion focused on the variability of synchrotron emission, 
the gamma-ray data is not that relevant, and has there¬ 
fore be only partially considered. Therefore, we fixed B 
to 0.038 G, derived from SED fitting of TeV gamma-ray 
data in Abdo et al. (2011). SED data were selected with 
the condition that XRT and optical/UV data were simul¬ 
taneously taken within a day. With this selection, we ob¬ 
tained 72 data points from our observations. To account 
for inverse Compton scattering, we used GeV gamma-ray 
data observed by the Fermi Gamma-ray Space Telescope. 
Abdo et al. (2011) reported that the GeV gamma-ray vari¬ 
ability is small, so we assumed that the gamma-ray flux 
was constant in our htting. In fact, GeV gamma-ray vari¬ 
ability coincident with the radio band has been reported 
(Lico et ah, 2014). In this study, however, gamma-ray 
variability was not essential for the fitting process, be¬ 
cause our fitting focused on the variability of synchrotron 
radiation. Therefore, we performed SED fitting under the 
assumption of a constant GeV flux. We used the 2FGL 
catalog value from 2008-2010 observations for our fitting 
(Nolan et ah, 2012). 

Figure 4 shows the light curves for X-ray and optical 
data. We performed SED fitting for different X-ray flux 
states, corresponding to the dates labeled in Fig 4. The 
one-zone SSC model provides a general description of the 
SED for both data in periods 1 and 2. However, upon 
closer examination, the one-zone SSC model may be in¬ 
sufficient for explaining simultaneously optical-UV and X- 
ray variability. Specifically, these results suggest that the 
optical-UV and X-ray emission may be coming from dif¬ 
ferent jet regions. We first discussed the simplest one- 
zone SSC model. Next, we will discuss the multi-emission 
model. The parameter values derived by SED fitting are 
similar to the values reported in Abdo et al. (2011). The 
variability in the Doppler factor and the co-moving blob 
radius R' were small (e.g., the variability of the Doppler 
factor was ^ ±20%). To explore the variability in these 
periods in detail, we also performed SED fitting using a 
fixed parameter for the Doppler factor D = 21 and co¬ 
moving blob radius R' = 5.2 x lO^^cm (Abdo et al., 2011), 
because these parameters did not show large variability 
in either of the two periods. Hereafter, we described the 
results of SED fitting with hxed parameters for B, D, and 
R'. Figure 5 and table 1 show the SED fitting results in 
periods 1 and 2. As a result, the SEDs well described the 
synchrotron component. Next, we discuss the variability 
mechanism of synchrotron emission based on these fitting 
results. In period 1, a SED with a higher synchrotron- 
peak frequency tended to have a higher luminosity. In 
contrast, in period 2, a SED with a lower synchrotron- 
peak frequency tended to have a higher luminosity. These 
results suggest that the variability mechanisms were dif¬ 
ferent for periods 1 and 2 . 

The variability of the synchrotron peak frequency to 
higher energy in period 1 was well described by the 
scenario of high-energy electron injection into the blob. 
Fitting results of the break energy of the electron distribu- 


Table 1. Fitting parameters with the fixed parameters 
B, D, and R' in period 1 and period 2 (Figure 5). 


Date 

pl“ 

p2^ 

Ke‘^ 

it 


MJD 55209 (A) 

2.39 

3.91 

103.4 X 10^*^ 

19.0 X 

T(F 

MJD 55230 (B) 

2.42 

3.87 

350.0 X lO'^o 

10.9 X 

10 ^ 

MJD 55252 (C) 

2.43 

4.80 

7.5 X 10^° 

50.8 X 

10 ^ 

MJD 55312 (D) 

2.46 

4.51 

5.0 X 10^° 

50.7 X 

10 ^ 

MJD 55537 (E) 

2.33 

4.54 

352.6 X 10^*^ 

12.2 X 

T(F 

MJD 55568 (F) 

2.34 

4.80 

500.0 X 10^° 

10.8 X 

10 ^ 

MJD 55598 (G) 

2.50 

4.56 

354.7 X 10^° 

9.5 X 

10 ^ 

MJD 55623 (H) 

2.37 

4.80 

500.0 X 10^° 

10.4 X 

10 ^ 


a: Index of the electron distribution (7 < 7(,), b: Index of the 
electron distribution (7 > 7 ;,), c: Normalization of the electron 
distribution, d: Break energy of the electron 


tion also indicate high-energy electron injection in period 
I. A scenario with energy injection by internal shocks in 
relativistic shells (shock-in-jet model) is one way to ex¬ 
plain such variability on a daily or monthly basis (e.g., 
Spada et al., 2001; Zhang et ah, 2002). The high peak 
in the synchrotron radiation could be obtained by the in¬ 
jection of energy into the emitting regions, and a con¬ 
sequent increase in the number of high-energy electrons. 
Polarimetric observation supports this scenario. The rel¬ 
atively high polarized emission in period I is well repre¬ 
sented by the shock-in-jet model, which naturally results 
in a high PD, due to blob compression. The shock-in¬ 
jet model also predicts the polarization vector orientation 
parallel with respect to the projected jet with compressed 
shocks. As mentioned previously, the position angle of the 
parsec-scale jet was measured to be —20 deg, a value close 
to the average PA of 150 deg in period I. The reason for 
the small variability of optical flux in period I compared 
with that of X-rays was that the optical emission comes 
from electrons with less energy than 7 ^. Decreasing the 
normalization value of the electron distribution suggests 
that the variability of the long-term component could be 
initiated by decreasing the number of low-energy elec¬ 
trons. 

The spectral variability in period 2 showed quite differ¬ 
ent behavior from that in period I. The variability in the 
X-ray band was relatively weak, compared with that of 
period 1. From the right panel of Figure 5 and Table 1, 
we were able to estimate the variability mechanism in pe¬ 
riod 2. One remarkable feature of the variability in period 
2 was the increase in the normalization of the electron 
distribution. This finding suggests that only the number 
of low-energy electrons increased, while changes in the 
spectral shape of the electron distribution were small in 
period 2. To reproduce such variability, we proposed a 
scenario in which insufficient acceleration or electron dif¬ 
fusion had a greater influence in period 2 , compared with 
that in period 1. Superposition of the multi-emission re¬ 
gion indicated weak polarized emission in period 2. In 
this case, electron diffusion was facilitated by weak con¬ 
finement. Ushio et al. (2010) discussed the possibility of a 
diffusive acceleration of electrons with subluminal shocks. 
A turbulent multi-emission zone model provides another 
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Fig. 4. Light curves for X-ray and optical flux. The dotted line indicates dates corresponding to SED fitting results shown in this pa- 
per. 




Fig. 5. SEDs of Mrk 421 with SSC model fitting with the fixed parameters of B, D and R' in pe¬ 
riod 1 and 2. Colors indicate the four periods. The parameter values are shown in table 1. 


possible explanation (e.g., Marscher, 2014). As mentioned 
above, SED fitting with the one-zone model is limited in 
its ability to simultaneously present optical-UV and X-ray 
variability. The multi-emission model also explains the 
small variations in polarization in period 1 (also see Lico 
et ah, 2014). In this case, the superposition of turbulent 
plasma downstream of the jet may explain the random 
behavior of the polarization. This model also explains the 
increase in the number of low-energy electrons, under the 
assumption of an increase in the diffuse emission region 
downstream of the jet. 

The SED model can also be constructed with a weak 
magnetic field (e.g., B < 0.038) instead of assuming a 
strong beaming factor. However, we cannot verify which 
model is more suitable, due to the coupling parameter. To 
distinguish between these models, temporally dense mon¬ 
itoring of this source is required for spectral evolutions 
predicted using a timescale of a few hours. Thus, tem¬ 


porally and energetically dense monitoring of this source 
would enable a more detailed study, using two-component 
models and time-dependent SSC scenarios. 

5. Summary and Conclusion 

To summarize our multi-wavelength observations of 
Mrk 421, systematic differences in the polarization and 
variability mechanisms were observed between periods 1 
and 2. The systematic component of polarization indi¬ 
cated the emergence of a long-term component in period 
1 with an aligned magnetic field. The variability of the 
whole-band spectra in period 1 was explained by the sce¬ 
nario of high-energy electron injection, such as that pro¬ 
posed by the shock-in-jet model, which also reproduced 
the polarized emission and direction of the PA. The vari¬ 
ability observed in period 2 can also be described by high- 
energy electron injection; however, in this case, the elec- 


























trons were not accelerated to ultra-relativistic velocities. 
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